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ABSTRACT 

In  an  earlier  paper  (Model  Atmospheres  (I)),  a  procedure  was  developed 
for  determining  the  most  probable  vertical  temperature  profile  associated 
with  the  occurrence  of  1^  global  temperature  extreme  at  mandatory-pressure 
levels  at  stations  in  the  North  American  Arctic.   The  same  technique,  based 
upon  a  variation  of  the  stepwise  multiple  regression  procedure,  was  employed 
in  the  current  study.   Whereas  the  radiosondes  investigated  in  Model  Atmos- 
pheres (I)  consisted  entirely  of  "checked-data"  quality,  those  stations 
designated  for  study  in  this  work  required  a  much  more  refined  data-screen, 
due  to  lack  of  initially  checked  radiosonde  report  quality.   Nevertheless, 
after  application  of  various  acceptability  criteria,  the  radiosondes  at 
each  station  were  arranged  in  the  same  format  as  employed  in  Model  Atmos- 
pheres (I).  There  remained  in  each  case  a  suitable  sample  population  to 
provide  significant  results.   The  ensuing  multiple  regression  analysis 
applied  to  the  geographically  and  climatologically  diverse  set  of  stations 
of  the  current  study  led  to  realistic  estimates  of  the  temperature  profiles 
which  were  conditionally  dependent  upon  the  existence  of  1^  extreme  forcing- 
level  temperature  T  at  previously  designated  pressure  levels  p  . 
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1.   Introduction 

The  present  work  frequently  to  be  abbreviated  in  what  follows  as  Model 
Atmospheres  (II)  is  a  sequel  to  a  research  report  Model  Atmospheres  (I)  by 
Martin  (1972) .   The  latter  report  dealt  with  the  development  of  vertically- 
consistent  temperature  profiles  which  are  contingent  upon  the  simultaneous 
existence  of  a  known  temperature  (at  the  1^  probability  level)  at  certain 
designated  stations  and  for  correspondingly  known  elevations »  The  complete 
set  of  stations  and  levels  for  which  these  tentative  MIL-STD-210  temperature 
extremes  were  considered  to  exist  are  shown  in  Table  1  after  Sissenwine 
(1970,  personal  communication)  for  both  the  cold-  and  warm-extreme  cases. 
The  extreme -sounding  envelopes  are  also  shown  in  Fig.  1  for  January  and  July 
climatology. 

Model  Atmospheres  (I)  dealt  with  those  extreme  sounding  sites  and  levels 
of  Table  I  indicated  by  asterisk  symbols,  which  stations  were  located  in  the 
North  American  Arctic  Basin,  Multiple  regression  procedures  were  developed 
for  specifying  the  vertical  temperature  profiles  at  these  North  American 
Arctic  stations.  All  of  the  stations  analyzed  in  Model  Atmospheres  (I) 
were  characterized  by  radiosonde  data  of  high  quality  and  by  a  type  of 
quality  control  termed  "checked  data"  processing.   Furthermore,  all  of  these 
Model  Atmosphere  (I)  cases  had  temperature-height  listings  in  the  concise 

pressure-level  format  p,  ,  k=l, ,21,  with  p,  ranging  from  1000  mb  to 

100  mb  as  shown  in  Table  2,  This  concise  arrangement  of  the  T(p) -profiles 
made  possible  a  simple  system  for  application  of  the  regression  procedures 
so  successfully  developed  in  this  earlier  report.  Model  Atmospheres  (I) 
thus  served  as  a  useful  pilot  project  for  the  development  of  mean  tempera- 
ture profiles  for  those  stations  not  indicated  by  the  asterisk  symbol  in 
Table  1,  which  station-data  has  become  the  basis  of  the  analysis  of  the 
present  paper. 
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Fig.  1.   Preliminary  version  of  the  MIL-STD-210B  Atmosphere  adopted  for 
use  in  this  study  (after  Sissenwine,  1970). 


TABLE  1.   Locations  of  proposed  extreme  temperatures  in  the 
cold  world-wide  and  warm  world-wide  cases  (after  Sissenwine, 
1970). 


Cold-Extreme 

Warm-Extreme 

Level 

Location 

January 
Mean 
(°C) 

STD. 
DEV. 
(°C) 

1^ 

Extreme 

(°C) 

Location 

July 
Mean 

STD. 
DEV. 

Extreme 
(°C) 

SFC 

Ojmjakon, 
USSR 

-50.6 

4.4 

-60.6 

Insalah, 
Algeria 

40.9 

3.9 

49.0 

850mb 

Ojmjakon, 
USSR 

-35.1 

5.5 

-47.0 

Insalah, 
Algeria 

28.9 

2.2 

34.0 

700mb 

Hall  Beach,* 
NWT 

-27.3 

6.5 

-42.4 

Baghdad, 
Iraq 

17.0 

2.2 

22.1 

500mb 

Resolute,* 

Nwr 

-43.0 

4.6 

-53.2 

New  Delhi, 
India 

-4.3 

3.9 

4.0 

300mb 

Thule,* 
Greenland 

-60.6 

2.2 

-66.0 

New  Delhi, 
India 

-25.8 

3,9 

-16.0 

200mb 

Thule,* 
Greenland 

-59.7 

7.2 

-75.0 

Alert,* 

NWT 

-42.3 

2.5 

-36.5 

150mb 

Karachi, 
Pakistan 

-65.0 

6.0 

-79.0 

Alert,* 

NWT 

-43.3 

2.5 

-37.5 

lOOmb 

Singapore 

-81.0 

2.2 

-86.0 

Thule,* 
Greenland 

-43.8 

2.2 

-37.2 

TABLE  2.  Arrangements  of  sounding  temperatures 
and  heights  by  pressure  levels  k=l,...,2,  in  each 
sounding. 


k  =  1 

2 

3 

4 

5 

k=6 

lOOOmb 

950 

900 

850 

800 

750 

k  =  7 

700 

650 

600 

550 

450 

450 

k  = 

12 

k  =  13 

400 

350 

300 

250 

200 

175 

k  = 

18 

k  =  19 

150 

125 

100 

k  =  19 

20 

21 

The  organization  of  Model  Atmospheres  (II)  proceeds  first  with  the  re- 
maining January  cold-extreme  cases  of  Table  1.   In  the  second  phase,  the 
July  warm-extreme  cases  are  considered.  No  attempt  is  made  in  this  paper 
to  restate  specific  results  of  Model  Atmospheres  (I) . 

Grateful  acknowledgements  to  the  Commander,  Naval  Weather  Service  Command 
should  be  made  at  the  outset  for  providing  all  of  the  data  tapes  for  the 
stations  listed  in  Table  1  from  the  original  data-tape  archives.  The  capable 
technical  assistance  of  the  Naval  Weather  Service  Environmental  Detachment, 
Asheville,  North  Carolina,  is  also  acknowledged.  All  data-samples  provided 
were  based  upon  the  period  1967-70  inclusive. 

Certain  problems  of  data  inaccuracy  occurred  in  connection  with  all  of 
the  soundings  of  the  present  study  in  both  the  summer  and  winter  cases.  The 
elimination  of  data  inaccuracies  is  the  subject  of  Section  2. 
2.   Data  processing. 

The  primary  data  deficiencies  in  the  soundings  considered  in  this  paper 
as  contrasted  with  those  of  Model  Atmospheres  (I)  were 

(i)   the  relative  scarcity  of  complete  soundings  particularly  due  to 
missing  data  below  the  850  mb  level,  and  above  the  200  mb  level. 

(ii)   inconsistency  in  the  number  of  reporting  levels: in  some  cases, 
complete  soundings  included  only  data  at  mandatory  levels.   In  others,  the 
number  of  temperature-reporting  levels  seemed  excessive  beyond  possible 
significance. 

(a)  Testing  and  development  of  vertically- cons is tent  soundings. 

In  maximizing  the  information  content  of  individual  soundings, 
invaluable  assistance  was  volunteered  by  Mr.  Donald  C.  Schertz  of  the  Naval 
Environmental  Prediction  Research  Facility  of  Monterey,  California. 
Mr.  Schertz  subjected  all  of  the  taped  soundings  to  a  computer-editing 


program  designed  to  scan  each  sounding  iteratively  between  reported  levels. 
The  main  objective  was  to  arrive  at  hydros tat ically  consistent  values  of 
the  reported  virtual  temperatures  and  heights  at  indicated  pressure  levels. 
In  doing  this,  superadiabatic  lapse  rates  above  850  mb  were  eliminated. 
Finally,  the  originally  listed  sounding  and  the  edited  sounding  were  printed 
out  in  succession  so  that  apparent  errors  in  the  original  were  easily 
recognizable.   Even  if  the  change  listed  by  the  editing  program  was  not 
accepted,  the  indicated  change  pointed  out  the  location  of  the  problem  data- 
level.   In  this  case,  maximum  weight  was  attached  to  the  hydrostatic  thick- 
ness equation    - 

Z   -  Z  =  -^  [1  -  0  ^  ^'^   ]  (1) 

where  p.,  p.,  are  successive  pressure  levels  in  the  vertical,  T  .  is  the 

A    B  VA 

virtual  temperature  at  the  base  of  the  layer  being  tested,  and  7    is  the 

lapse  rate  of  virtual  temperature  between  p.  and  p  »  The  virtual  temperature 

A       B 

T       at  the  upper  level  is  then  given  by 
VB 

The  conversion  between  T(p)  and  T  (p)  is  given  by 

/l_ 


T^(P) 

1  -H   -^^^^^^   exp[21.656   • 

5418.1 

where  T  (p)  is  the  listed  dewpoint  temperature  (  K)  at  level  p.   The  denom- 
inator of  (3)  reduces  essentially  to  unity  when  p  <  500  mb. 

When  T(p)  was  available  at  p  =  200  mb  but  not  at  higher  levels,  a 
"bootstrap"  regression  procedure  was  utilized  to  generate  T(IOO).   This 
was  done  by  utilizing  a  regression  equation  for  T(IOO)  in  terms  of  predictors 


at  all  uniform  temperature- levels  (similar  to  those  of  Table  2)  below  p  = 
200  mb.   This  procedure  generally  led  to  good  specification  of  T(IOO)  for 
the  dependent  sample,  and  also  to  reasonable  values  of  T(IOO)  when  the 
latter  was  originally  listed  as  missing. 

For  any  sounding  case  not  discarded  by  gross-error  checks,  if  a  tempera- 
ture T(p)  was  missing  at  a  designated  level  p,  of  Table  2,  the  temperature 
T(p  )  was  interpolated  using  the  hydrostatically-consistent  equations  (1), 

(2),  (3)  with  p  >  p  >  p„.   This  procedure  was  also  applied  to  the  T(p,  ) 
A    k    B  K 

values  between  200  and  100  mb,  when  the  latter  had  been  generated  by  the 
regression  method.   In  this  case,  it  was  also  necessary  to  develop  a  value 
Z(IOO)  first  in  order  to  apply  Eq.  1.  This  was  done  by  choosing  as  a  first 
guess  for  the  lapse  rate  the  station  sample-mean  height  Z  (100),  where 

^  T(200)  -  T(IOO) 

^'^      Z   -  Z 

^100   200 

While  Eq.  4  introduces  some  error  in  the  lapse  rate  in  the  layer  (100,200), 
the  thickness -computation  Az (100,200)  tends  to  be  error-compensated.  The 
resulting  height  computations  are  then  used  only  in  computation  of  T(p,  ) 
values  at  Z(p,  )  by  Eq.  2  to  fill  out  the  temperature  matrix  of  Table  2. 
(b)  The  case  of  missing  surface  reports. 

For  the  most  part,  the  major  effort  of  this  study  hinges  on  the 
use  of  surface  temperature  reports  to  determine  the  vertical  structure  above 
the  surface.   In  approximately  25fo   of  all  cases  listed  in  the  station  files, 
the  surface  temperature  report  was  missing.   In  nearly  all  these  cases, 
the  entire  sounding  was  discarded,  since  the  temperature  structure  between 
850  mb  and  the  surface  was  usually  extremely  variable.  Two  exceptions 
occurred. 


(i)  New  Delhi  January  soundings.   Here  the  temperature  profile  to  be 
determined  was  based  upon  a  known  initial  value  T  (150).  The  regression- 
generated  profile  could  be  determined  in  a  large  nuiriber  of  cases  down  to 
850  mb.  For  the  smaller  subset  of  sounding  cases  containing  surface  reports, 
only  surface  means  and  standard  deviations  were  computed  since  these  mean 
results  indicated  the  absence  of  extreme  temperatures  at  the  surface  level. 

(ii)  Ojmjakon,  USSR,  January  soundings.  For  this  station,  the  sur- 
face inversion  was  so  strong  as  to  be  a  relatively  consistent  feature  of 
the  vertical  profile  below  850  mb.   Since  the  1000  mb  heights  were  almost 
always  listed,  it  was  considered  useful  to  determine  the  1000  -  850  mb 
thickness,  in  the  form 


T  (1000,850)  =  Z(850)  -  Z(IOOO)  ^3^^ 


^d  /,    1000  \ 
rV""    "850  ) 


and  then  to  equate  the  result  of  (5a)  to  the  height -weighted  mean  tempera- 
ture defined  by      - 

_    [Z(850)  -  726.](^^^^|^^t_Is)  +  [726,-  Z(IOOO)]  T^ 
^  "  Z(850)  -  Z(IOOO)  ^^^^ 

For  Ojmjakon,  the  surface  elevation  Z  =  726.0  meters  was  employed,  and 

the  terrain  contribution  to  the  mean  temperature  T  was  made  isothermally 

at  the  value  detennined  from  Eq.  (5b)  for  T 

The  values  of  T  resulting  from  Eq.  5  showed  good  agreement  with  the 

listed  surface  mean  and  standard-deviations  of  Table  1.  The  computation 

of  T  was  highly  dependent  upon  Az (850, 1000)  and  this  parameter  had  to  be 
s 

carefully  checked  for  invalid  values.  Such  values  usually  arose  from  re- 
ported errors  in  Z(IOOO),  which  could  be  detected  from  an  analysis  of  the 
time-sequencing  of  the  Z(IOOO)  reports,  (i.e.,  by  gross-error  checks  of 

Z(IOOO)). 
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After  eliminating  such  thickness  errors,  the  T  -values  were  accurate 
enough  to  be  specified  at  a  highly  significant  level  in  terms  of  the  over- 
lying temperature  variables  in  the  format  of  Table  2. 

As  noted  previously,  virtually  all  other  stations,  e.g.^  Insalah  and 
Baghdad,  required  definition  of  the  surface  temperature  to  consider  the 
possibility  of  surface  extremes.  Hence,  if  the  surface  report  were  missing, 
the  entire  sounding  was  discarded.  At  Singapore,  the  data-quality  was 
excellent  and  no  data  rejections  were  necessary. 

The  use  of  New  Delhi  January  data  was  attempted  solely  as  a  substitute 
for  Karachi,  where  preliminary  indications  (Table  1)  led  one  to  expect  a 
cold  extreme  at  150  mb. 

As  a  result  of  the  data-processing  outlined  here,  the  following  station 
sounding  sets  arranged  in  the  format  of  Table  2  were  developed: 

Winter  Sounding  data. 
Ojmjakon  Singapore  New  Delhi  (for  Karachi) 

169  cases  86  cases  158  cases,  p  <  850  mb, 

92  cases,  surface  and  above 

Summer  Sounding  data. 
Insalah  Baghdad  New  Delhi 

72  cases  92  cases  159  cases 

3.   Development  of  regression-generated  profiles,  January  extremes 

The  stepwise  regression  procedure  of  BIMED  02R  was  utilized  in  conjunc- 
tion with  the  checked  data  files  for  each  winter  case  station  of  this  study. 
The  procedure  is  virtually  identical  to  that  described  in  Model  Atmospheres 
(I).  The  objective  at  each  station  was  to  develop  by  a  modified  least- 
squares  technique  the  temperature  specification  for  each  T   (temperature  at 
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a  mandatory  level),  where  the  temperature  T  successively  assxjmes  each  of 
the  indexed  variables  (from  Table  2)  listed  below,  apart  from  the  forcing- 
level  temperature. 

T- ,   the  temperature  at  the  surface 

T, ,   the  temperature  at  850  mb 

Ty,   the  temperature  at  700  mb 

T^^,   the  temperature  at  500  mb 

T-c>   the  temperature  at  300  mb 

T,,,   the  temperature  at  200  mb 

T^Q,   the  temperature  at  150  mb 

T  -,   the  temperature  at  100  mb 

The  procedure  then  was  to  develop  a  least-squares  fit  for  the  specifi- 
cation of  T,,(P^,)  in  the  form 
M  M 

where  T  (p  )  is  the  temperature  at  the  pressure  level  p  where  the  extreme 
temperature  listed  in  Table  1  occurs.  The  variable  X^  represents  a  linear 
combination  of  up  to  7  independent  temperature  variables,  not  including  the 
forcing-level  temperature  T  .  The  variable  X^  involves  temperature-predictors 
which  are  highly  correlated  with  the  dependent  variable  Z,.  The  variable  T 
is  forced  into  (6)  at  each  step  of  the  eight-step  regression  equation  sought. 

An  analysis  in  the  form  of  Eq.  6  will  be  applied  to  each  of  the  follow- 
ing January  station  data  files: 

(A)  Ojmjakon,    Tj  =  Tj(Sfc) 

(B)  Ojmjakon,    Tj  =  Tj(850) 

(C)  Singapore,   T^  =  Tj(lOO) 

(D)  Singapore,   T^  =  Tj(150) 

(E)  New  Delhi,   Tj  =  Tj(150) 
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The  regression-generated  statistics  for  cases  (A) , (B) , , . . , (E) ,  include 
means,  standard  deviations,  multiple  correlation  coefficients  R  and  standard 
errors  of  estimate.  The  primary  test  of  the  specification  procedure  is 
through  the  least  squares  foirmula  (6)  which  in  general  yields  an  R  close  to 
unity  and  a  standard  error  much  reduced  from  the  standard  deviation.  Three 
analyses  have  been  conducted  in  each  case  and  these  analyses  are  summarized 
in  tabular  form  for  each  case. 

These  analyses  are  based  essentially  upon  the  regression  procedure  applied 
to  (a)  the  full-sample  of  T  -data,  (b)  the  same  regression  to  the  10^  nominal 
extreme  sample,  and  (c)  to  the  1^  subsample  extreme  of  T  -values.  A  nominal 
extreme  may  be  defined  for  each  station  and  T,  according  to  the  Gaussian 
probability  percentiles 

^.Ol^Pj*  =  ^J  ±  2.3267  aj 

where  the  minus  sign  applies  for  the  cold  extreme.  Values  of  T  ,a  are  taken 

from  the  full-sample  analyses  in  part  (a)  of  each  table.   In  each  of  tables 

a,b  is  shown  also  the  regression  statistics  R  and  a  ,  where 

£ 

R  is  the  multiple  correlation  coefficient 


with  a  the  standard  error,  after  application  of  the  specification  equation 
for  T^. 

It  is  to  be  noted  that  in  tables  3, 4,..., 7  below  that  there  is  little 
systematic  deviation  in  the  standard  errors  resulting  from  treatment  (b) 
relative  to  that  for  the  full-sample  treatment  (a).  These  results  give 
credibility  to  the  requirement  of  a  meaningful  forcing-level  temperature  T^. 
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However,  in  most  cases,  the  requirement  of  a  1^  extreme  of  T  does  not  yield 
a  significantly  different  vertical  profile,  level  for  level,  from  that  re- 


quired by  the  existence  of  a  10%  extreme  of  T  .  This  appears  to  be  due  to 
the  surprisingly  small  standard  deviations  at  TL^-levels  for  the  stations 
under  consideration  here,  which  is  in  sharp  contrast  to  the  rather  large 
standard  deviations  in  the  Model  Atmospheres  (I)  study, 
3(A).  Ojmjakon,  T  (Sfc). 

The  results  for  the  three-subs ample  stratifications  are  shown  in 
Tables  3a,b,c.   It  is  to  be  noted  that  the  1^  extreme  listed  here  is  -61.3  C 
which  is  colder  than  that  of  Table  1.   However,  in  part  (b) ,  the  20fo   sub- 
sample  has  been  analyzed  rather  than  the  10^  subsample  normally  to  be  treated 
under  part  (b).   This  is  due  to  the  emphasis  (AFCRL,  1973)  placed  upon  the 
"coldest  feasible  operating"  surface  temperature  of  -51  C,  which  corresponds 

to  T  «-(Sfc)  at  Ojmjakon.  A  consistent  deviation  between  the  T,,(c)  and  T,.(b) 
.20    /     J  J  MM 

profiles  does  appear  in  Table  3  in  that  the  low-level  temperatures  are 
slightly  colder  in  the  former  case,  and  become  somewhat  warmer  at  high 
elevations  (p  <   150  mb).  The  essence  of  these  tabular  comparisons  may  be 
visualized  by  examining  Fig.  2. 
3(B).   Ojmjakon,  Tj(850. 

A  similar  set  of  comparisons,  level  for  level,  may  be  made  by  examina- 
tion of  Table  4  (a,b,c).  The  graphs  of  Fig.  3  show  that  the  T  Q^i^^O)   atmo- 
sphere tends  to  remain  colder  than  the  mean  at  all  levels  below  p  =  150  mb, 
but  then  tends  to  become  slightly  warmer  at  p  =  100  mb. 
3(C).   Singapore,  Tj(lOO). 

Here  we  examine  the  winter  extreme  atmosphere  models  from  the  tropo- 
spheric  top  rather  than  from  the  base.   It  is  to  be  noted  that  observed  ifo 
extreme  T  (Table  5)  is  colder  than  that  listed  in  Table  1.  Fig.  4  illustrates 
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Fig.    2.      Regression-generated  mean  January  temperature  profiles  over  OJMJAKON, 
USSR,   corresponding  to  the  forcing-level  temperature  T    (Sfc),      The  heavy  solid 
profile  is  based  on  the  ifo  cold  extreme  sample  at  the  surfaceo      Both  profiles 
are  based  upon  values   listed  in  Table  3. 
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Fig.  3.  Regression-generated  mean  January  temperature  profiles  over 
OJMJAKON,  USSR,  corresponding  to  a  forcing-level  temperature  T  (850). 
The  heavy  solid  profile  is  based  on  the  1^  cold  extreme  sample  at  850mb. 
Both  profiles  are  based  upon  values  listed  in  Table  4o 
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Fig.  4,      Regression-generated  mean  January  temperature  profiles  over 
Singapore  corresponding  to  a  forcing-level  temperature  T   (100).      The 
heavy  solid  profile  is  based  on  the   1^  cold  extreme  sample  at   100  mb. 
Both  profiles   are  based  on  values   listed  in  Table  5. 
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the  type  of  modification  which  may  be  expected  between  the  T  (c)-  and 

T», (a) -profiles,  with  the  latter  based  upon  the  full-data  sample.  The 

M 

former  case  is  colder  than  the  mean  profile  at  Singapore  to  about  300  mb, 
and  below  this  level  becomes  somewhat  warmer. 
3(D).   Singapore,  Tj(150). 

Table  6  is  relevant.  The  test  conducted  here  was  done  primarily 
because  the  Karachi  data  was  unavailable  in  any  useful  voltime  or  quality. 
The  sample  mean  [case  (a)]  is  cold  for  this  level  in  agreement  with 
meridional  profiles  of  the  climatological  data  (Crutcher,  1970)  in  the 
region  between  the  equator  to  25N  in  January.  The  assignment  of  the  1^ 
cold  extreme  at  150  mb  listed  in  Table  6  is  characteristic  of  a  cooler 
than  normal  upper  troposphere,  but  case  (c)  shows  a  tendency  toward  a 
warmer  mid-troposphere  (e.g.,  at  500  and  700  mb).  The  nature  of  any  pro- 
file deviations  closer  to  ground  level  is  not  clear  from  the  statistics 
of  Table  6.   Fig.  5  shows  the  comparison  between  the  mean  and  "extreme" 
atmospheres.  Another  interesting  result  is  that  the  case  (c)  extreme  at 
150  mb,  T  (150)  =  -74.3  °C  corresponds  to  T  (100)  =  -82.0  °C  which  is 

considerably  warmer  than  the  climatological  extreme  T  _. (100)  =  -87.1  C 

•  u  J. 

noted  in  Section  3(D). 
3(E).   New  Delhi,  T  (150). 

Again  this  case  was  tested  primarily  as  a  substitute  for  T  ^-(150), 
originally  listed  in  Table  1  as  occurring  at  Karachi.  However,  there  was 
no  reliable  data  at  Karachi.   The  part  (a,b,c)  regression  results  for 

New  Delhi  are  listed  in  Table  7  for  100  mb, ,850  mb,  levels  only.  The 

listed  I5S  extreme  T  ^^^(150)  =  -70.7  °C  at  New  Delhi  is  far  warmer  than 
that  of  Table  2  for  Karachi  (-79  C) .  The  latter  extreme  presumably  has 
been  incorrectly  estimated  because  of  the  large  standard  error  attributed 
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Fig,  5.  Regression-generated  mean  January  temperature  profiles  over 
Singapore  corresponding  to  a  forcing-level  temperature  T  (150),  The 
heavy  solid  profile  is  based  on  the  1^  cold  extreme  sample  at  150  mb, 
Both  profiles  are  based  on  values  listed  in  Table  6, 
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Fig,  6.  Regression-generated  mean  January  temperature  profiles  over  New 
Delhi  corresponding  to  a  forcing-level  temperature  T  (150) .  The  heavy 
solid  profile  is  based  on  the  1^  cold  extreme  sample  at  150  mb.   Both 
profiles  are  based  on  Table  7.  The  extension  below  850  mb  is  included 
in  Table  7. 
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to  Karachi  at  150  mb,  which  seems  to  be  due  to  unreliable  instrxxmentation. 
The  1^  extreme  Tj(150)  =  -70.7  °C  observed  here  for  New  Delhi  (Table  7) 
is  considerably  warmer  than  that  found  at  Singapore,  and  therefore  Singapore 
should  qualify  in  place  of  Karachi  as  approximately  representative  of  a 
world-wide  cold  extreme  at  150  mb. 

The  surface  regression  at  New  Delhi  has  not  been  carried  out.   It  is 
clear  from  the  listed  means  and  standard  deviations  at  the  surface  (Table  7) 
that  the  diurnal  variation  tends  to  obscure  interdiumal  variations.  This 
presumably  is  a  characteristic  of  the  winter  monsoon  regime  at  this  loca- 
tion. Fig.  6  indicates  that  the  vertical  profile  for  a  New  Delhi  cold 
extreme  T  q. (150)  is  associated  with  a  cool  upper  troposphere  (relative  to 
the  mean)  with  a  reversal  toward  a  warm  lower  troposphere  below  300  mb. 
3(F).   Summary  of  cold  extreme  regions 

Recalling  the  results  of  Model  Atmospheres  (I)  in  conjunction  with 
the  present  study,  we  can  draw  two  well-defined  distinctions.  The  coldest 
temperatures  at  100  mb  and  possibly  at  150  mb  are  to  be  found  in  equatorial 
regions.  The  coldest  temperatures  between  700  mb  and  200  mb  are  found  in 
extreme  conditions  in  the  Arctic  circle  under  cold  vortex  (aloft)  condi- 
tions, this  temperature  regime  being  somewhat  more  extreme  than  that  found 
at  Ojmjakon.  The  latter  station,  however,  is  consistently  cold  aloft  at 
all  of  these  levels,  whereas  the  Arctic  circle  stations  tend  to  undergo 
periodic  "explosive  warmings"  even  in  January  above  500  mb.  Of  course,  the 
coldest  surface  stations  (up  to  850  mb)  are  found  in  the  landlocked  areas 
of  Siberia. 
4.  Development  of  regression-generated  warm  profiles,  July  extremes. 

The  regression  procedure  of  Section  3  was  followed  here  with  the  July 
data-files  of  Insalah,  Baghdad  and  New  Delhi  using  appropriate  forcing-level 
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temperatures  T  in  deriving  specification  equations  for  each  mandatory  level 
T  y  T  o   In  general  a  discussion  of  the  results  in  the  order  T  (Sfc),  T  (850), 
....,  T  (300)  will  be  given  below. 
4(A).  Insalah,  Tj(Sfc). 

The  results  of  the  regression  procedure  on  the  full-sample  data  and 
the  nominal  10^  extreme  sample  are  given  in  Table  8.   In  this  case,  T  q, (Sfc) 
was  not  observed;  in  fact  was  not  really  closely  approached.  This  was  due 
to  the  fact  that  the  sotinding  record  of  the  first  level  retrieved  averaged 
approximately  40  meters  above  the  station  elevation.  At  this  level,  the 
amplitude  of  the  diurnal  temperature  was  much  reduced  so  as  to  be  unrepre- 
sentative of  the  Table  1  listing  of  the  surface  warm  extreme  temperature. 
The  T  _^ (Sfc)  =  43.7  C  value  listed  under  Table  8(c)  is  associated  with  a 
1^  histogram  count  from  the  actual  distribution  of  "surface"  temperatures. 
Fig,  7  shows  that  the  T»t(c)  mean  profile  is  considerably  cooler  in  the  upper 

troposphere  than  is  T^,(a). 

M 

4(B).  Insalah,  Tj(850). 

Here  the  suggested  extreme  of  Table  1  was  observed,  as  is  evident 

from  T^(c)  =  34.7  °C  of  Table  9(c).  The  T  (c)-profile  which  was  correlated 

with  this  T  ^. (850)  then  turns  out  to  be  consistently  warmer  than  both  T„(b) 
.Ui  M 

and  T^(a),  level  for  level.  The  fact  that  T  (c)  is  warmer  than  T  (a)  is 
clearly  depicted  in  Fig,  8.  The  most  unusual  feature  of  the  t(c) -profile 
is  the  relatively  cool  value  at  the  surface  T^(Sfc)  =  41.7  C,  indicating  a 
possible  association  with  a  local  storm  effect. 
4(C).  Baghdad,  Tj(Sfc) . 

Although  Baghdad  is  suggested  as  having  the  extreme  warm  value  at 
700  mb  (Table  1) ,  it  turned  out  that  the  Baghdad  data  listed  surface 
temperatures  as  hot  as  that  suggested  for  the  1^  warm  extreme  (Table  1) 
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occurring  in  the  Algerian  desert.  Thus  the  data  for  Baghdad  was  also  used 
for  Tj.(Sfc)  since  that  for  Insalah  was  inconclusive  in  this  respect.   It 
turned  out  that  the  radiosonde  surface  temperatures  were  well  coordinated 
with  shelter  temperature  values,  and  two  reports  of  49°C  did  occur  in  this 
area  of  the  Iraqi  desert.  The  corresponding  profiles  are  listed  in  Table 
10  for  the  (a),  (b),  (c)  stratifications.  Fig.  9  is  indicative  of  the  varia- 
tion between  the  extreme  and  the  full-san?)le  case,  and  shows  only  slightly 
warmer  temperatures  in  the  former  case  at  all  mandatory  levels.  As  noted 
in  some  previous  comparisons  the  difference  between  the  10^  and  1%   "extreme" 
cases  seems  to  be  insignificant. 
4(D).   Baghdad,  Tj(700),  0000  GMT 

Here  the  suggested  nominal  value  of  T  _.  (700)  given  by  Eq.  7  was 
exceeded  upon  two  separate  occasions.  Oddly,  this  occurred  only  with  the 
0000  GMT  reports,  the  warmest  1200  GMT  T  (700)  reports  being  somewhat  cooler. 
Hence,  a  separate  analysis  of  the  two  different  sets  of  sounding  times  at 
Baghdad  was  made  and  is  summarized  in  Tables  11,  12.  The  tendency  for  a 
warm  extreme  T  „^ (700)  to  occur  at  night  suggests  that  nighttime  subsidence 
may  be  occurring  at  the  700  mb  level  following  the  maxiratim  insolation  during 
the  preceding  day.  Fig.  10  shows  that  the  tropopause  temperature  tends  to 
be  cooler  aloft  under  these  same  conditions,  i.e.  T    (100)  <  t^(lOO)  <  T^  (100), 
and  that  this  same  trend  is  in  effect  all  all  upper  levels  p  <   200  mb.  Thus 
following  a  day  of  insolational  maximum,  the  statistics  suggest  the  concomitant 
occurrences  of  high-level  tropospheric  convergence  with  700-mb  subsidence. 
Both  of  these  features  tend  to  support  the  differences  between  profiles  (a), 
(b)  and  (c). 
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Fig,  7.   Regression-generated  mean  July  temperature  profiles  over  Insalah, 
Algeria,  corresponding  to  a  forcing-level  temperature  T  (Sfc),  The  heavy 
solid  profile  is  based  on  the  1^  extreme  warm  sample  at  the  surface.   Both 
profiles  are  based  on  Table  8. 
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Fig.  8.   Regression-generated  mean  July  temperature  profiles  over  Insalah 
corresponding  to  a  forcing-level  temperature  T  (850).  The  heavy  solid 
profile  is  based  on  the  1^  extreme  warm  sample  at  850  mb.  Both  profiles 
are  based  on  Table  9. 
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Fig.  9.   Regression-generated  mean  July  temperature  profiles  over  Baghdad, 
Iraq,  corresponding  to  a  forcing-level  temperature  T  (Sfc)  at  1200  GMT. 
The  heavy  solid  line  is  based  on  the  1^  extreme  warm  sample  at  the  surface. 
Both  profiles  are  based  on  Table  10. 
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4(E).   Baghdad,  T  (700),  1200  GMT. 

The  results  in  this  case  are  listed  in  Table  12,  with  an  extreme  warm 
T  (700)  =  19.7  °C,  which  is  lower  than  the  0000  G3^  extreme.  However, 
the  case  (c)  surface  temperature  T^  (Sfc)  =  49.0  C  corresponds  to  the 
daytime  T  (Sfc)  extreme.  Table  12  and  Fig,  11  both  show  the  T^(c)  pro- 
file to  be  slightly  warmer  at  all  levels  than  the  \r(b)  and  T^(a)profiles 
at  1200  GMT.   In  contrast  with  the  0000  GMT  profile-extreme,  the  1200  GMT 
is  dominated  by  convection  at  all  tropospheric  levels  in  view  of  the 
strong  surface  layer  lapse.   In  fact  this  case  resembles  that  of  T  (850) 
at  Insalah,  apart  from  the  surface-level  perturbation  of  temperature  at 
Insalah. 

The  Baghdad  1200  GMT  data  file  revealed  no  occurrence  of  a  warm 
extreme  at  850  mb  which  equalled  that  observed  over  the  Algerian  desert 
(Insalah).  Consequently,  the  problem  of  generating  a  model  atmosphere 
T^(c)  based  upon  the  contingent  condition  of  a  T  ^^ (850)  at  Baghdad  was 
not  considered.   In  fact  comparison  of  Tables  8  and  10  show  that  apart 
from  the  surface  and  700  mb  levels  the  Insalah  midday  profile  at  all  other 
levels  is  warmer  than  that  over  Baghdad.  The  apparent  difference  at  the 
surface  has  been  explained  in  Section  4(A),  and  the  explanation  of  a 
700-mb  extreme  at  Baghdad  has  evidently  to  do  with  stronger  anticyc Ionic 
control  aloft  which  manifests  itself  clearly  in  the  0000  GMT  soundings. 
4(F).  New  Delhi,  Tj(500)  and  Tj(300). 

These  cases  should  first  be  reviewed  separately  as  summarized  in 
Tables  13  and  14,  respectively.   It  should  be  noted  that  the  proposed 
T  -(500)  =  4  °C  of  Table  1  was  not  observed  at  New  Delhi.   This  Table  1 
overestimate  seems  to  be  due  to  an  overestimate  of  the  standard  deviation 
in  conjunction  with  Eq.  7.  However,  it  is  significant  to  note  that  Table  13 
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Fig,  10,  Regression-generated  mean  July  temperature  profiles  over  Baghdad 
corresponding  to  a  forcing-level  temperature  T  (700)  at  0000  GMT.  The 
heavy  solid  profile  is  based  on  a  1^  extreme  warm  sample  at  700  mb.  Both 
profiles  are  based  on  Table  11, 


36 


60 

to 

ij 

C 

u   .u 

0) 

•1-1 

ra    u 

CO 

w 

•a    to 

3 

1      (X 

t-l 

g 

«k 

!—(        .• 

CO 

3  ^-v 

^ 

«4-l  O 

o 

0) 

V 

o  r>. 

B 

rt 

jz  ^-^ 

<u 

V4 

iJ  H 

u 

M 

4J 

O  *« 

X 

• 

«^ 

■U    O 

o 

^-v 

TS 

o 

rt 

W    <U'fe'l> 

r^ 

•X3 

u  •-* 

1— ( 

1 

x 

Q  a. 

r>~ 

60 

»+-i    E 

a; 

VO 

<U     CO 

V4   en 

1-t 

4J 

^-N    <u 

o 

to 

CO 

CO  e 

4J 

v_x    (y 

13 

o 

to 

0) 

W 

60 

1-4 

*J     4J 

P 

<y 

>J      X 

•r-l 

> 

r:    o 

•O 

•H 

o 

^  e 

C 

1—1 

j_i 

O 

^ 

CO 

a. 

5 

c; 

•    > 

CO 

}-l 

t-i 

0) 

c; 

3 

u 

V4 

C_( 

w 

s^^*?. 

u 

S 

n 

*J  o 

o 

u 

0) 

CO   ,-i 

o 

^ 

u 

o 

a 

O    r-4 

w 

o 

ex.    CO 

a> 

CM 

>» 

E    C 

^-1 

i-l 

u 

O    v-l 

3 

o 

AJ     H 

iJ 

E 

4J 

o 

CO 

o 

C3 

c 

^ 

M 

73 

<y 

IW 

a 

'D 

cu 

r3 

>.  J2 

f1 

G) 

e 

r-t      u 

o 

^1 

2 

4-1 

CO 

•W 

-)   o 

ro 

u 

c 

0) 

T-* 

CO 

w 

O     03 

OJ 

(U 

O 

>  u 

B 

•H 

O     (U 

•T3 

u 

w 

•r-l 

.u 

^    «4-l 

rH 

w 

1     <u 

(U 

•A 

^^ 

> 

a) 

o 

•H  ^-x 

.— 1 

,— ( 

u 

o  ^ 

1 

w 

V^  ^~' 

60 

to 

o 

c 

4J 

c 

«+-!     AJ 

•H 

^ 

o 

^1 

•o 

o 

T^ 

<D     CO 

c 

Pi, 

w 

X    CU 

3 

(U 

w 

iJ 

O 

V* 

<y 

•  m, 

tO- 

;4 

VI     o 

i-l 

60 

CO  •-• 

(U 

i-< 

<1) 

D.  jr 

<: 

PJ 

^  e 

iJ 

O    CO 

• 

• 

o   w 

TO 

^-^ 

CN 

r^ 

iJ 

o 

r-) 

s^ 

CO 

o 

H 

-r^ 

r^ 

W 

.— t 

v^ 

(-3 

a    >. 

H 

^ 

t-i 

/"-N 

•"^S 

o 

u^ 

H 

H    ^ 

V-' 

O 

«w  u 

oo 

• 

•       ^ 

-o 

W    AJ 

u 

>-l     CO 

CO 

w  w 

CO 

a 

• 

m 

• 

•  *w 

CD 

u 

U    <44 

u 

.— * 

U     Qi 

3 

o  o 

<!■ 

S 

CJ  o 

a 

^ 

• 

• 

X! 

"O 

>  o 

^— ' 

AJ 

oo 

CO 

Q 

TO 

to 
CO 
o 

o 
11 


H 
Al 


CO 

so 


oo 

•  «^ 
U  U 
U    to 

w  w 


• 

• 

«M 

4J 

u 

m 

.-1 

u 

4) 

3 

o 

o 

X 

o 

o 

t3    >  O 
CO  Q 


CO 

o  o 

XO 


> 

o 


c^ 


CN 


cn       t-i 


CM 


CM 

I 


CM 


m 


00 


St         VO         o 

.^^       i/^       r^ 

I  I         I 


VO 
00 


m 


<t 

m 

P^ 

^ 

ir> 

<t 

VO 

CO 

r^ 

CJ\ 

r-4 

o 

I-l 

CO 

VO 

o 

CO 

t» 

<r 

<TN 

<T> 

O* 

C3N 


r-l 

VO 

r-l 

o 

to 

cy\ 

CX) 

■CM 

a^ 

CTi 

o> 

a\ 

CT> 

C7^ 

ty» 

a\ 

CX) 

CO 

00 

ON 

CM 

r- 

T-I 

<3N 

o 

VO 

r-l 

<t 

r^ 

m 

VO 

VO 

vO 

CO 

CO 

-d- 

VO 

CO 

I-l 

a\ 

• 

• 

• 

• 

• 

• 

• 

• 

<SI 

<N 

o 

m 

I-l 

r-l 

CM 

CM 

r- 

CM 

CM 

o 

1^ 

(Nl 

r^ 

CM 

r^ 

r- 

vO 

r^ 

in 

VO 

to 

r* 

• 

• 

• 

• 

♦ 

• 

• 

• 

VO 

o 

CO 

CO 

vO 

VJ- 

VO 

r-l 

•<1- 

CO 

r-l 

1 

CN 
1 

1 

lO 
1 

1 

CO 

CM 

CJv 

<!■ 

VO 

CO 

to. 

CM 

r^ 

CO 

ON 

I-l 

<y\ 

cy> 

r~ 

CM 

o 

<t 

I-l 

r-l 

liO 

<t 

^ 

Si; 

CO 

ON 

CM 

-:!• 

to 

CM 

ON 

•sJ- 

r>^ 

CM 

vO 

<J\ 

vf 

r^ 

c^ 

ON 

00 

00 

C3N 

CJN 

ON 

o\ 

ON 

ON 

r-l 

o 

r-4 

00 

CO 

o 

CO 

cr. 

CO 

ON 

VO 

(y\ 

00 

VO 

C^N 

CO 

O 

<t 

-* 

00 

r-l 

I-^ 

o 

lO 

CO 


CM 


CM 


CM 


CM 


CM 


CM 


00 

lO 

00 

00 

CO 

CM 

VO 

00 

r-l 

o 

ON 

r-l 

1^ 

<!• 

CO 

00 

• 

• 

• 

• 

• 

• 

• 

• 

CO 

r- 

CO 

CM 

VO 

to 

00 

CM 

<t 

CM 

r-l 

• 

CM 

1 

1 

in 
1 

1 

^ 

S 

^ 

m 

to 

O 

o 

O 

o 

o 

o 

o 

ON 

lO 

o 

o 

o 

o 

to 

o 

<JN 

00 

r^ 

to 

CO 

CM 

t-H 

r-4 

37 


20  t- 


-  200 


-  300 


-  500 


-  700 


-  850 


1000 


Fig.  11,  Regression-generated  mean  July  temperature  profiles  over  Baghdad 
corresponding  to  a  forcing-level  temperature  T  (700)  at  1200  GMT.   The 
heavy  solid  profile  is  based  on  a  1^  extreme  warm  sample  at  700  mb  at 
1200  GMT.   Both  profiles  are  based  on  Table  12. 
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Fig.  12.  Regression-generated  mean  July  profiles  over  New  Delhi  corres- 
ponding to  a  forcing-level  temperature  T  (500) .  The  heavy  solid  profile 
is  based  upon  a  1^  extreme  warm  sample  at  500  mb.   Both  profiles  are  based 
on  Table  13 « 
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Fig.  13.  Regression-generated  mean  July  temperature  profiles  over  New 
Delhi,  corresponding  to  a  forcing-level  temperature  T  (300).  The  heavy 
solid  profile  is  based  on  the  1^  extreme  warm  sample  at  300  mb.   Both 
profiles  are  based  on  Table  14. 
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does  indicate  a  shift  in  warmest  July  500  nib  temperatures  from  Baghdad  to 
New  Delhi.   This  phenomenon  of  T  qj^(500)  at  New  Delhi  is  evidently  associated 
with  the  summer  monsoon  circulation  over  New  Delhi  and  its  associated  release 
of  latent  heat  at  the  forcing  level  T  (500), 

J 

In  the  case  of  an  extreme  forcing-level  temperature  at  the  300  mb  level, 
latent  heat  release  goes  on  to  much  higher  levels.  Hence  the  T  (b),  T^(c) 
profiles  of  Table  14  are  warmer  above  300  mb  than  those  of  Table  13.  These 
features  of  the  "extreme"  New  Delhi  soundings  are  credible  in  view  of  the 
stramer  monsoon  effect.  Note  that  in  both  extreme  profiles  of  Fig.  12  and 
Fig.  13,  respectively,  the  surface  level  "extreme"  data  are  cooler  than 
those  of  the  full-sample  means,  presumably  because  of  reduced  surface  pre- 
cipitation in  the  full-sample  case  (a). 
5.  Conclusions 

The  data  treated  in  this  study  is  in  general  based  upon  less  sounding 
data  than  that  of  Model  Atmospheres  (I).  As  a  result  in  addition  to  check- 
ing the  data  for  hydrostatic  consistency,  they  have  been  checked  for  also 
meteorological  reality  wherever  possible.  In  general,  the  locations  of 
extreme  temperature  sites  are  located  as  originally  suggested  (Sissenwine, 
1970,  personal  communication).  More  specifically,  some  of  the  original 
ifo   extremes  of  Table  1  appear  to  have  been  overestimated  (warm  cases)  or 
underestimated,  depending  upon  the  validity  of  the  standard  deviation. 
The  data-processing  employed  in  this  study  seems  to  have  generated  more 
realistic  standard  deviations  as  well  as  means. 

A  further  effort  along  the  lines  of  Model  Atmospheres  (I)  and  the 
current  study  is  anticipated  for  a  set  of  Oceanic  Station  Vessels  and 
for  several  island  stations  where  oceanic  extremes  (rather  than  world-wide) 
are  known  to  occur. 
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